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We study theoretically the spatial correlations between the intensities measured at the input 
and output planes of a disordered scattering medium. We show that at large optical thicknesses, 
a long-range spatial correlation persists and takes negative values. For small optical thicknesses, 
short-range and long-range correlations coexist, with relative weights that depend on the optical 
thickness. These results may have direct implications for the control of wave transmission through 
complex media by wavefront shaping, thus finding applications in sensing, imaging and information 
transfer. 
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I. INTRODUCTION 


The study of wave scattering in disordered media is 
an active field of research, stimulated both by innovative 
applications in imaging and sensing [l| and by fundamen¬ 
tal questions in mesoscopic physics [ 3 - In the last few 
years, the possibility to control the propagation of op¬ 
tical waves in complex media in the multiple scattering 
regime has been demonstrated using wavefront shaping 
techniques Si- This breakthrough offers new perspec¬ 
tives for imaging and communication through complex 
media im. The initial schemes make use of optimiza¬ 
tion techniques requiring intensity measurements in the 
transmitted speckle, which in terms of practical applica¬ 
tions is a serious drawback. Finding a way to control 
the transmission and focusing of light through a strongly 
scattering medium from measurements of the reflected 
speckle only is an issue of tremendous importance. Pro¬ 
gresses have been made recently by taking advantage of 
the memory effect [I®, with imaging capabilities lim¬ 
ited to relatively small optical thicknesses. Nevertheless, 
the connection between the reflected and the transmit¬ 
ted speckle patterns generated by a disordered medium 
in the multiple scattering regime has not been addressed 
theoretically so far. 

In this paper, we make a step in this direction by study¬ 
ing theoretically and numerically the statistical correla¬ 
tion between the intensities measured in the transmitted 
and the reflected speckle patterns. The spatial intensity 
correlation function C{y,y') is defined as 


C{y,y') 


{6I{r)SI{Y')) 
(J(r)) (/(rO) 


( 1 ) 


where the notation (...) denotes a statistical average over 
disorder and SI{y) = /(r) — (/(r)) is the intensity fluc¬ 
tuation. This correlation function has been extensively 
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studied in the context of wave scattering and mesoscopic 
physics PH. Theoretical approaches often make use 
of the canonical slab or waveguide geometries (for a re¬ 
view see P, P, Pi and references therein), where either 
transmitted or reflected intensity is considered, or con¬ 
sider point sources in an infinite or open medium and 
compute intensity correlations for two points inside or 
outside the medium PP|. It seems that the intensity 
correlation function for two points lying on different sides 
of a slab medium has not been studied, and that the exis¬ 
tence of a correlation has only been mentionned through 
passing PP. In this work, we study the correlation 
between the intensities in the input and output planes 
of a strongly scattering slab, as sketched in Fig. [T] Us¬ 
ing numerical simulations and analytical calculations, we 
show that for optically thick slabs a correlation persists 
and takes negative values. Moreover, at smaller opti¬ 
cal thicknesses, short and long-range correlations coexist, 
with relative weights that depend on the optical thick¬ 
ness. We believe these results to be a step forward for 
the control of transmission through strongly scattering 
media, thus finding applications in sensing, imaging and 
information transfer. 

The spatial intensity correlation in a speckle pattern 
can be split into three contributions, historically denoted 
by Cl, C2 and C3 P,P|, as follows: 

C(r, r') = Cl (r, r') + C2 (r, r') + C3 (r, r'). ( 2 ) 

The first term Ci corresponds to the Gaussian statistics 
approximation for the field amplitude, and is a short- 
range contribution, whose width determines the aver¬ 
age size of a speckle spot P. C2 and C3 are non- 
Gaussian long-r ang e correlations that decay on much 
larger scales H, In the diffusive regime, and for two 
observation points lying on the same side of the scattering 
medium, these three contributions have different weights 
such that Cl C2 C3. In the reflection/transmission 
configuration considered here, and at large optical thick¬ 
ness, we show that an intensity correlation persists and 
is dominated by C2 because of the short-range behavior 
of Cl. In this case, the long-range character of C2 is con- 
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fered by its algebraic decay with respect to the distance 
between the two observation points. Moreover, this cor¬ 
relation is negative, a result that may have implications 
in the context of wave control by wavefront shaping. For 
smaller optical thicknesses, we also show that a crossover 
can be found between regimes dominated by Ci and C 2 , 
respectively. 

The paper is organized as follows. In Section m we 
study the reflection/transmission correlation function at 
large optical thicknesses. First, we use numerical simu¬ 
lations to compute the correlation function without ap¬ 
proximation and to describe its main features. Second, 
we present the analytical calculation of the C 2 contribu¬ 
tion to the reflection/transmission correlation function in 
the multiple scattering regime, for both two-dimensional 
and three-dimensional geometries, and show that C 2 is 
the leading contribution at large optical thicknesses. In 
Section Ell we study the correlation function at small 
optical thicknesses, a regime in which the Ci contribu¬ 
tion dominates. In Section lYl as a consequence of the 
negative value of the correlation that is found in the mul¬ 
tiple scattering regime, we describe some peculiarities 
of the statistical distribution of reflected/transmitted in¬ 
tensities. Finally, in Section |3 we summarize the main 
results and discuss some implications for the control of 
wave transmission through disordered media. 


II. REFLECTION/TRANSMISSION 
CORRELATION AT LARGE OPTICAL 
THICKNESS 

A. Numerical analysis 

In this section we present exact numerical simulations 
of wave scattering in the multiple scattering regime. We 
restrict ourselves to a 2D geometry for the sake of com¬ 
puter memory and time. We consider a slab of scattering 
material, characterized by its thickness L and its trans¬ 
verse size D (we keep D > 6L in order to avoid finite size 
effects), as depicted in Fig. [TJ Our purpose is the calcu¬ 
lation of the correlation function C{rfi,rT), where is 
a point located on the input surface (reflection) and tt 
is a point located on the output surface (transmission). 

1. Method 

To proceed, we use the coupled dipoles method [ 2 ^ 
to calculate numerically the intensity in the transmitted 
and reflected speckle patterns. Repeating the calcula¬ 
tions for a large number of conhgurations of disorder (po¬ 
sitions of scatterers) allows us to compute statistics. The 
system contains N randomly distributed non-overlapping 
point scatterers, and is illuminated by a plane wave from 
the left at normal incidence. We deal with TE-polarized 
waves with an electric field oriented along the invariance 
axis of the system (scalar waves). The resonant point 



FIG. 1: (Color online) Schematic representation of the scat¬ 
tering medium. The slab of thickness L and transverse size D 
is illuminated from the left by a monochromatic plane-wave at 
normal incidence. The correlation between the reflected and 
transmitted speckle patterns is characterized by the correla¬ 
tion function between the intensities at points r_B (reflection) 
and Tt (transmission). 

scatterers are described by their electric polarizability 

2r 

= ” 12 ?- 

k^lyui — Wo -I- *F/2) 

where wq is the resonance frequency and F the linewidth. 
This specific form of the polarizability fulfils the optical 
theorem (i.e. energy conservation). From the polarizabil¬ 
ity the scattering cross section cr(w) = fco|Q;(w)p/4 and 
the scattering mean-free path £(w) = [ptT(w)]“^ can be 
deduced, where p = N/{LD) is the number density of 
scatterers. In the following, we consider scatterers at 
resonance (w = wq) in order to reach large optical thick¬ 
nesses with a reasonable number of scatterers (typically 
a few hundreds). In the coupled dipoles formalism, the 
exciting field Ej on scatterer number j is written as 

N 

Ej = Eo {vj) + a (lu) kl ^ Go (r^ -Yk)Ek (4) 

where Gq is the 2D free-space Green function given by 
Go(r — r') = (i/4) Hg^^(fco|r — r^|), being the Han- 
kel function of first kind and order zero. Equation (|4|) 
dehnes a set of N linear equations that are solved by a 
standard matrix inversion procedure. Once the exciting 
field is known on each scatterer, the field E{r) and the 
intensity /(r) = |E(r)p at any position r inside or out¬ 
side the scattering medium can be calculated by a direct 
summation, using 

N 

E{r) = Eo (r) -b a (w) ^ Gq (r - r^) Ej. (5) 
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2. Numerical experiment 

We have carried out numerical simulations in the mul¬ 
tiple scattering regime with an optical thickness b = 
Lji = 7. This choice of optical thickness is limited by 
the number of configurations that can be calculated in 
a reasonable computer time in order to get a sufficiently 
accurate statistics to compute averaged values (requiring 
typically 10® configurations). The correlation function 
C'num obtained from numerical calculation is plotted in 
Fig. [5] (red solid line) versus the lateral shift Ar between 
the observation points in the reflected and transmitted 
speckles (see the geometry in Fig.[T]). 



FIG. 2: (Color online) Reflection/transmission intensity cor¬ 
relation C'num given by the numerical simulations (red solid 
line) and analytical correlation C 2 given by Eq. (1281) (blue 
dashed line). Multiple scattering regime with b = 7 and 
koi = 10. 

Surprisingly an intensity correlation subsists even for 
large optical thicknesses (multiple scattering regime). 
Moreover the reflection/transmission correlation function 
at large optical thickness takes a negative value around 
Ar = 0. This means that {SI{rii)SI{rT)) < 0, showing 
that the probability to have a dark spot in the trans¬ 
mitted speckle in lateral coincidence with a bright spot 
in the reflection speckle (and vice-versa) should be sub¬ 
stantial. This property, that might have implications for 
the control of wave transmission by wavefront shaping, 
is investigated more precisely in section [TV] 


B. Analytical derivation of the C 2 
reflection/transmission correlation function 

To get more insight on the numerical result presented 
above, we present the calculation of the C 2 contribution 
to the intensity correlation function for scalar waves in 
both two-dimensional (2D) and three-dimensional (3D) 
geometries. Intuitively, we expected C 2 to be the lead¬ 
ing contribution at large optical thicknesses, since C 3 is 
always negligible compared to C 2 and Ci vanishes ex¬ 


ponentially with the optical thickness in this particular 
reflection/transmission configuration. 

For a dilute system such that koi ^ 1, where k = 
ujjc = 27r/A with c the speed of light and A the wave¬ 
length in vacuum, the intensity correlation functions 
can be calculated analytically using a perturbative ap¬ 
proach. Since the calculation for reflection-reflection 
or transmission-transmission correlations is detailed in 
textbooks or review articles we do not give 

all details here but rather focus on the specificity of 
the reflection-transmission geometry. For the analytical 
derivation, we consider that the transverse size D of the 
slab is infinite. 

1 . Average intensity 

To compute the correlation function C 2 , we first need 
to compute the average intensity. The starting point is 
the Bethe-Salpeter equation that reads [1, [l3| 

{E{r,u;)E*{r',uj)) = {E{r,uj)) {E*{r',u;)) 

+ J (G(r,ri,a;)) (G'*(r',r 2 ,w)) 

X r(ri, r2, ra, r4, uj) (£'(r3, uj)E*{r 4 ,u})) dridr2dr3dr4. 

( 6 ) 

In this equation, (G'(r, r',w)) is the average Green 
function that links the average field {E{r,ui)) to a 
source dipole p located at position r' via {E{r,uj)) = 
{G{r,r',uj)) p. The operator r(ri, r 2 , ra, r 4 , w) is the 
irreducible vertex that contains all multiple scattering se¬ 
quences connecting (ra, r 4 ) to (ri, r 2 ). The exact cal¬ 
culation of this complex object is out of reach, but an 
approximate expression to first-order in the small param¬ 
eter l/{koi) and for independent scattering, known as the 
ladder approximation, can be derived M- Omitting the 
frequency dependence for simplicity in the following, it 
reads 

r(ri, r2, ra, r4) = 75(ri - r2)(5(r3 - r4)^(ri - ra) (7) 
where 

4fco/.^ for 2D TE waves (2D), 

Att/£ for 3D scalar waves (3D). 

Plugging Eq. ([7]) into Eq. ([HI) yields 

(/(r)) = |(E^ + 7 J\{G{ry))\^{I{r'))dr'. (9) 


Physically, this closed equation means that only contri¬ 
butions for which E and E* follow the same scattering 
path have a significant weight in the average intensity (all 
cross terms vanish). The first term /s in Eq. ([H]) is the 
ballistic (or coherent) component of the average inten¬ 
sity and the second term /d is the diffuse part. In order 
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to get explicit expressions for these two quantities, we 
need to compute the average Green function, or equiva¬ 
lently the average field. As a consequence of the Dyson 
equation [l^, in the independent scattering limit, the av¬ 
erage Green function obeys a propagation equation in an 
effective homogeneous medium, defined by an effective 
wavevector fceff. We have 


(G'(r)) - < exp(ifceff|r|) 

[ 47r|r| 

and {E{r)) = Eq exp(zfceffz) 


(2D), 

(3D), 


( 10 ) 

( 11 ) 


where k^g = ko + to first order in l/(fco^) and z 

is the coordinate along the direction normal to the slab. 
The ballistic intensity is readily deduced: 


Ib{z) = /oexp[-z/£]. (12) 

Regarding the diffuse intensity Id{y)^ we can rewrite 
Eq. ([9]) in the following way 


/D(r) = 7 j\{G{v,v'))\^ [Ib{v') + Id{v')\Av' . (13) 


An analytical expression of the diffuse intensity can be 
obtained in the diffusive limit where |r — r^| '> Making 
use of the translational invariance of the medium, the 
Fourier transform of Eq. (na) reads 


where 


-iD(q) = 7^(q) [-^^(q) + -fr’(q)] 


A{q) = - X 


1 I y/T+<^ 


7 I arctan(g£) 

I 


(2D), 

(3D). 


(14) 


(15) 


As we consider large distances, we can perform a second 
order Taylor expansion of 1/A(q) for qi 1 which leads 
in the real space to a diffusion-type equation 

-^Aloir) = iBir) (16) 

where d € {2, 3} is the space dimension. Solving Eq. (flHl) 
in a slab geometry, we obtain 


Ind 


+ {z + zo) (l - y) exp 


{L + zo — z) ^1 -I- 


L 


-/odexp(^-|^ (17) 


where zq is the extrapolation length needed to account 
for the boundary conditions at the input and exit surfaces 
of the slab [2^, For a dilute index-matched slab, its 
expression is given by [13, 


7r£/4 (2D), 

2£/3 (3D). 


This expression of the diffuse intensity Ib is a priori valid 
for distances z such that z ^ £, for which the diffusion 
approximation holds, but it surprisingly gives reasonably 
reliable results even for z <£. Adding the ballistic term, 
it also gives reliable results for the full average intensity 

ij{z)) = 1 b{z) A Id{z) (19) 

even for relatively small optical thicknesses. 


2. Long-range C 2 contribution 


The intensity correlation function is a fourth order cor¬ 
relation in terms of field amplitude. Physically, a correla¬ 
tion is created when the two pairs of fields that constitute 
the intensities in the correlation function share a com¬ 
mon history in the scattering process. Regarding C 2 , the 
crossing occurs during the propagation of the intensities 
inside the system and is described by a complex object 
known as a Hikami vertex , and denoted by H in the 
following. The propagation of the intensity between the 
slab surfaces and the crossing is described by the ladder 
operator, denoted by L in the following. The expression 
of C 2 is given by [ij, [13, [H, Ull 

C2ivR,YT) = -j-r -,—/ dridr2dr3dr4 

{I{yr)) (/(rr))] J 

X dpidp2dp3dp4|(G(rfl, rz))H(G(rT, Yi))\^ 

X L{^2, P2)LiT^i, Pi)H{pi, P2, P3, P4) 

xT(p3,r4)L(p3,r3)K£;(ri))n(E(r3))P (20) 

which can be rewritten diagrammatically in the following 
form: 


C2{rR,rT) = 


rR 


( 21 ) 


Tt 



where the hatched box is the Hikami vertex and the other 
boxes are the ladder operators. Thick lines represent the 
average Green function and thick dashed lines stand for 
the average electric field. 

The ladder operator is defined as 

L(r,r') = 7'^(r-r')+7 J \{G{r,r"))\'^L{r",r')dr'' (22) 
and is represented diagrammatically as 


O OK) OK>K3 


L(r,r') = 


O OK) OKKKD 


+ ... (23) 


Zo = 


(18) 
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where circles and thick horizontal solid lines represent 
scattering events and average Green functions, respec¬ 
tively, the top line standing for the field amplitude E 
and the bottom line for its complex conjugate E*. Thin 
vertical solid lines link scattering events involving iden¬ 
tical scatterers. 

An analytical expression of the ladder operator can 
be obtained in the diffusive limit previously used for the 
computation of the average intensity. Making use of the 
translational invariance along the direction of the slab 
interface, the Fourier transform of the ladder operator 
with respect to transverse variables can be obtained from 
Eq. (1^ . and reads 


L(K,z,z) = — 

sinh[A:( 2 :< 2 : 0 )] sinh[A:(L + Zp - z>)\ 

sinh[Ar(T -|- 2 ^ 0 )] 

where z< = min(z, z'), z> = max(z, z'). 

The Hikami box can also be calculated in the limit 
» 1, and its expression reduces to [ll|, MM 

H{Pi,P2.P3,P4) = h I Y[S{p-p,)Vp^-Vp^dp (25) 

i—1 

where 

^ fC/(16.S) (2D), 

\C/(24»tS) (3D), 

To have an explicit expression of (72, we first compute 
the integrals involving ri and ra using 

J i(r,r')l(-6^(i‘'))l^dr'= 7 (/(r)) . (27) 

As we deal with large optical thicknesses, we can replace 
the average intensity by its diffuse component. Regard¬ 
ing the integrals involving r 2 and r 4 , we assume that the 
ladder operators vary slowly at the scale of the scatter¬ 
ing mean-free path i. This amounts to replacing r 2 by 
rji and r 4 by ry in the ladder positions. We end up 
with an explicit expression of the C 2 contribution to the 
correlation function, given by 


C2{rR,rT) = 



where 


{I{rn)) (/(rr)) Jo \ L J 

^ + l| sinh() 3 ') - gcosh(g) 

X sinh(gzo/L)^ 

{L + 2zoy sinh[g(l + 2zo/L)]^ ^ 


F 



16 cos{qAr/L) 



Trko£ q 
27 

p—Jo(gAr/L) 

^ /hQ-CXy 


(2D), 

(3D). 

(29) 


The expression of the C 2 contribution to the reflec¬ 
tion/transmission correlations function, together with its 
comparison to full numerical simulations, is the main re¬ 
sult of this paper. The dependence of the amplitude on 
the system size L is approximatively with d the di¬ 

mension of space. This shows that even at large optical 
thicknesses, a correlation subsists between the intensities 
measured in the transmitted and the reflected speckles, 
and that this correlation is dominated by a contribution 
of the C 2 type. 

It is important to keep in mind that this expression 
has been derived in the framework of the diffusion ap¬ 
proximation (in particular bulk average Green functions 
have been used). Its validity in the geometry consid¬ 
ered here where both reflected and transmitted intensi¬ 
ties contribute is checked by comparison to full numerical 
simulations in Fig. O It is interesting to note that the 
analytical and numerical calculations are in very good 
quantitative agreement, showing that the diffusion ap¬ 
proximation, used to derive the analytical results, is very 
accurate even for an optical thickness 6=7 (that is not 
very large) and a geometry involving a reflected intensity 
(that always involves short scattering paths). 

We have seen that C 2 decays algebraically with the 
system size L. Thus it is interesting to analyse the be¬ 
havior of the full correlation function at smaller optical 
thicknesses, and in the crossover between the multiple 
and single scattering regimes. This problem can be ad¬ 
dressed numerically, by using the numerical method de¬ 
scribed previously in various scattering regimes, and is 
the subject of the following section. 


III. REFLECTION/TRANSMISSION 
CORRELATION AT SMALL OPTICAL 
THICKNESS 


A. Numerical analysis 


We consider a small optical thickness 6 = 0.5, cor¬ 
responding to the single-scattering regime. The reflec¬ 
tion/transmission intensity correlation function calcu¬ 
lated numerically in this regime is shown in Fig. [3] (red 
solid line). The large oscillations observed on a scale on 
the order of the wavelength are not described by the C 2 
contribution, and are expected to be a signature of the 
short-range Ci contribution. To check the validity of this 
assumption, we have to compute the Ci contribution to 
the correlation function. At small optical thickness, a 
quantitative calculation would require to go beyond the 
diffusion approximation and to account properly for the 
boundary conditions [1^. Since our purpose in this sec¬ 
tion is only to support qualitatively the analysis of the 
general trends observed in the numerical simulations us¬ 
ing a simple model, we keep using the diffusion approxi¬ 
mation to estimate the Ci contribution. 
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FIG. 3: (Color online) Reflection/transmission intensity cor¬ 
relation Cnum given by the numerical simulations (red solid 
line) and analytical correlation including all terms C'i-I-C'i-I-C '2 
(blue dashed line). Single scattering regime with b = 0.5 and 
koi = 10 . 


Equation (1321) shows that Ci is the square of the scat¬ 
tered field correlation function. Indeed, Ci can also be 
seen as the correlation that woidd be observed for a 
field with Gaussian statistics, for which this factorization 
holds dH. Starting from Eqs. (l32l) and (l33l) . the usual 
way to derive the analytical expression of Ci consists in 
replacing r 2 by tt (for a speckle computed in transmis¬ 
sion) or r 2 by r/j (for a speckle computed in reflection) 
in the ladder positions. In the reflection/transmission ge¬ 
ometry, this simplification cannot be performed because 
the relative distance between the points and can 
be very large compared to the scattering mean-free path 
£. However, in order to get an explicit expression, we can 
make use of Eq. (EH) which leads to 

B{rR,rT)= f 7(G(r/i,r2)) (G*(rT,r2)) (/(r2))dr2. 

(34) 


B. Cl contribution and specificity of the 
refiection/transmission configuration 


The integration over r 2 is then performed using the 
residue theorem (the details are given in appendix ET). 
This leads to 


In the scattering sequences picture, the Gi correlation 
is created by interchanging the amplitudes between two 
independent ladders at the last scattering event. The Gi 
contribution reads 


CiirR^rr) = 777—77^777—77 [ dridr 2 dr 3 dr 4 

(J(rfl)) (/(rr)) J 

X (G*(rH,r4)) (G(rT,r4)) (G(r/i,r2)) (G*(rT,r2)) 

X T(r2,ri)L(r4,r3)|(E(ri))n(£;(r3))|^ (30) 

In terms of diagrams, Eq. dsni) can be rewritten as 



In Eq. (1^ . the scattering of both pairs of fields from 
points ri to r 2 and from r 3 to r 4 , respectively, is de¬ 
scribed by a ladder propagator. The mixing of am¬ 
plitudes at the last scattering event is represented by 
four different average Green functions. The integrals in 
Eq. (1501) can be factorized, leading to 


where 


Ci{rR,rT) 


\B{rR,rT)\‘^ 
(/(rfl)) (/(rr)) 


(32) 


B{yr,Yt) = j {G{tr,T2)) {G*{yt,V2)) 

X L(r 2 ,ri)|(£'(ri))pdridr 2 . (33) 


B{vR,rT) = ih 


exp [—{ik' + k'')L — iK. ■ Ar] 


4(fc'2 -H fc" 2 ) 
d{(l + zo/t)Mi + (I - zq/^) exp(-L/£)M2} 


L -(- 2 zo 
- (d - 1)M3 


dK 


. (27r) 


d-1 


(35) 


where Mi, M 2 , M 3 , k' and k” are given by Eqs. (TMI) . 
(lAQI) and (lAlOl) and Eq. (IA4I) . respectively. The last 
integral over K is performed numerically. 

In the specific geometry considered here, another con¬ 
tribution has to be added, in which the ballistic intensity 
contributes as one of the intensities involved in the corre¬ 
lation function [32|. Indeed, the ladder operator involves 
at least one scattering event, and does not account for sit¬ 
uations in which there is no scattering event before the 
field interchange. Such contributions to the correlation 
function can be important for small optical thicknesses, 
where the ballistic contribution is not negligible. This 
leads to a correction to the Gi correlation function, that 
we denote by G^, and whose expression is 


X J {G*{rR,r2)) (G(rT,r2))T(r2,ri)|(E;(ri))|^dridr2 

-|- c.c. (36) 


In terms of diagrams, the above expression can be rewrit- 
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ten as 



(37) 

Making use of the quantity tt) defined in Eq. (1551) . 
we obtain 


C[(rR,rT) 


2^[{E*{rR)) {E{rT)) B{rn,rT)] 
{I{rn)) (/(rr)) 


(38) 


It is important to note that the C'l contribution intro¬ 
duced here, and that involves the average field, is not 
negligible compared to the usual Ci contribution. This 
correction should be added to the Ci contribution when 
the product of the ballistic fields {E*{rii)) {E^tt)) can¬ 
not be neglected, as in the reflection/transmission geome¬ 
try at low optical thickness, or in the reflection/reflection 
geometry at any optical thickness. We can show from 
Eq. (1551) that the Ci and C[ contributions to the cor¬ 
relation function decrease exponentially with the optical 
thickness h = Lj 1. This behavior explains why C 2 dom¬ 
inates at large optical thickness. But in the single scat¬ 
tering regime, an important contribution of the Ci -I- C[ 
term is observed. This is clearly seen in Fig. [3] (blue 
dashed line), in which the sum of the three contributions 
Cl + C'l + C 2 calculated analytically in the same two- 
dimensional geometry used for the numerical simulation 
is plotted. Qualitatively, the behavior observed in the nu¬ 
merical simulation is fairly reproduced by the analytical 
approach. 

Moreover, since the Ci + C( contribution decays ex¬ 
ponentially with the medium thickness T, a crossover is 
expected towards a regime dominated by C 2 when the 
optical thickness increases. This also shows that the re¬ 
flection/transmission geometry studied here may be rel¬ 
evant to put forward experimentally the influence of the 
C 2 contribution (in the pure reflection or transmission 
geometries, the Ci contribution is always the leading con¬ 
tribution) . 


IV. STATISTICAL DISTRIBUTION OF 
REFLECTED AND TRANSMITTED 
INTENSITIES 

Analytical and numerical results at large optical thick¬ 
ness have shown that a correlation between reflected and 
transmitted intensities exists. Surprisingly, this correla¬ 
tion function takes negative values around Ar = 0. Hav¬ 
ing {SI{rR)SI{rT)) < 0 for Ar = 0 qualitatively suggests 
a high probability to have a bright (dark) spot in the 
transmitted speckle in coincidence with a dark (bright) 
spot in the reflected speckle. In order to address this 


question in more quantitative terms, we have studied the 
full statistical distribution of the intensities. More pre¬ 
cisely, from the numerical simulations, we have extracted 
the statistical distributions of the product of the fluctu¬ 
ating part of the intensities at Ar = 0, defined as 

SI{rR)6I{rT) 1 , , 

(/(r«))(/(rT))J ^ ^ 

whose average value is the intensity correlation C{rR, tr) 
at Ar = 0. 


P^SIrSIt) = P 



FIG. 4: (Color online) Statistical distribution of the 

product of the normalized reflected/transmistted intensities 
P^SIrSIt) at Ar = 0. Multiple scattering regime with 6 = 7 
and kol = 10. 


The statistical distribution P{5IrSIt) obtained for 
6 = 7 and koi = 10 is shown in Fig. [4] The distribu¬ 
tion exhibits an asymmetric shape, with a most likely 
value at SI{rR)SI{rR) = 0. Due to this asymmetric 
shape, the distribution cannot be simply characterized 
by its first moments. In particular, for 6 = 7 and 
koi = 10 , we find that the probability to have trans¬ 
mitted and reflected intensity fluctuations with opposite 
signs P{SI{rR)SI{rT) < 0) = 0.47, while we could have 
expected a much larger value (above 0.5) from a naive 
argument based on the negative sign of the correlation 
C{rR,rT)- Nevertheless, it is interesting to study statis¬ 
tics under some constraints. In particular, we have stud¬ 
ied the probability of having SI{rT) < 0 (a spot in the 
transmitted speckle darker than the average intensity) 
under the assumption that SI{rR) > n(/(r^)), with 
n G {0,1,2} (i.e. for a coinciding spot in the reflected 
speckle with increasing brightness). The results are sum¬ 
marized in Tabled] Interestingly, we find that the prob¬ 
ability P{6I{rT) < 0 ) increases substantially with the 
brightness of the reflected speckle spot, an information 
that is not contained in the intensity correlation function. 
In consequence, if the reflected intensity in a speckle spot 
is large compared to the average reflected intensity, the 
transmitted intensity in the coinciding spot in the trans¬ 
mitted speckle is smaller than the average transmitted in¬ 
tensity with a large probability. This result may have im¬ 
plications in the context of light focusing through opaque 











scattering media by wavefront shaping. Indeed, maximiz¬ 
ing the intensity in a reflected speckle spot might, with a 
high probability, leads to a minimization of the intensity 
in the corresponding transmitted spot. A precise study 
with optimized wavefronts (beyond plane-wave illumina¬ 
tion) is left for future work. 


b 

P 

Q 

q 

5I{rR) > 0 

SlivR) > {I{rn)) 

Slim) > 2(/(rfi)) 

0.5 

0.48 

0.57 

0.57 

0.62 

0.67 

1 

0.50 

0.61 

0.64 

0.69 

0.74 

2 

0.47 

0.64 

0.64 

0.67 

0.69 

4 

0.47 

0.65 

0.65 

0.67 

0.67 

7 

0.47 

0.66 

0.66 

0.67 

0.67 


TABLE I; Probability of having (i7(ri{)5/(rT) < 0 (p), of hav¬ 
ing SI{rT) < 0 (q) at different optical thicknesses (6) and un¬ 
der some constraints (5I(r_R) > n {I{ru)) with n € {0,1, 2}) 
for Ar = 0 and koi = 10. 


V. CONCLUSION 

In summary, we have studied analytically and numeri¬ 
cally the spatial correlation between intensities measured 
in the reflected and transmitted speckles generated by a 
slab of disordered scattering medium. We have demon¬ 
strated the existence of a reflection/transmission correla¬ 
tion. At large optical thicknesses, the spatial correlation 
persists and is dominated by the C 2 contribution, thus 
exibiting a long-range behavior. Interestingly, this corre¬ 
lation takes negative values. At small optical thicknesses, 
the correlation is dominated by a Ci-type contribution, 
which contains the usual Ci term and an additionnal 
term C[ involving the ballistic intensity. 

The statistical connection between transmitted and re¬ 
flected speckles might be of interest for wavefront shap¬ 
ing methods used to focus and image through scattering 
media. Since for practical implementations only the re¬ 
flected speckle can be measured and controlled, a knowl¬ 
edge of the probability to get a bright (dark) spot in the 
transmitted speckle in coincidence with a dark (bright) 
spot in the optimized reflected speckle could be a great 
advantage. As a first step towards this goal, we have 
studied the statistical distributions of reflected and trans¬ 
mitted intensities, and have identified situations in which 
the probability of coincidence of bright and dark spots on 
opposite sides of the medium is high. 

Finally, a refinement of the analytical model would be 
beneficial to deal with optical thicknesses for which the 
diffusion approximation fails to give quantitative results. 
One possibility could be to developp a semi-analytical 
approach (coupling analytical expressions and numeri¬ 
cal calculations) based on the radiative transfer equa¬ 
tion [ 3 ^ that described accuratly short and long scatter¬ 
ing paths [13 • 
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Appendix A: Analytical calculation of B 

In the reflexion/transmission configuration the usual 
approximation used to calculate Ci breaks. Indeed in 
transmission/transmission or in reflexion/reflexion one 
usually manages to separate in Eq. (l33l) integrals over 
r 2 and ri because of the small distance between the two 
points where we compute the correlation. In our config¬ 
uration, because of the large distance between these two 
points, we have to calculate explicitly 

Using the symmetry of the problem, taking the Fourier 
transform of the two Green functions over their trans¬ 
verse coordinates, we simplify Eq. (l33l) the following way: 


BivR, rr) = 1 7 {Iiz 2 )) (G(z 2 , K)) {G*{z 2 - L, K)) 

X exp[-fK • Ar]dz 2 

where Z 2 is the depth at which the two fields separate. 
To calculate the expression of the Green function in this 
mixed domain, one can see that: 


{G{z2,K)) = J (G(fc,K))exp[ffcz2]^ 

exp[ikz2] dk 

fc2+if2_fc2^27r' ^ ^ 

Equation (IA2I) can be calculated using the residue the¬ 
orem using the fact that k'^g = k^ + iko/i and noting 
that: 

1 _ 1 

“ {k-k+){k + k+) 

where 

fc+ = ^Jkl + iko/e- = k' + ik" 

is the pole with a positive imaginary part (i.e. k" > 0). 
Thus Eq. (IA2I) can be calculated with the residue theo¬ 
rem by integrating over a contour composed of the real 
axis and a upper-half-circle the radius of which tends to 
infinity. We obtain that 


(A3) 

(A4) 



(G(2:2,K)) = ^exp(fA:+2:2)- (A5) 
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Following the same steps we can show that With these expressions, we can rewrite the final ex¬ 
pression of rx) as a Fourier transform given by 

{G*{z2 - L, K)) = ^ exp[zfc;(z 2 - L)] (A 6 ) Eq. (051). 

leading to 


B{rR,rT) = 1 j {I{Z2)) 


exp[2iA:'2:2] 


4(^/2 fc // 2 ) 

X exp[—iK ■ Ar]d22 


exp[—ik’^L] 

dK 


(27r)" 


(A7) 


Using the expression of the average intensity we have to 
perform three integrations denoted respectively by Mi, 
M 2 and M 3 : 


Ml = / {L + Zq — Z 2 ) exp\2ik'Z2]dz2 

Jo 

1 -I- 2ik'{L + zq) — exp[2z/c'L](l + 2ik'zo) 
^ 4 p 2 ’ 

M2= {zq + Z 2 ) exp[2zfc'z2]d-22 

Jo 

2ik'zo -1-1- exp[2iA:'L][l — 2ik'{L + zq)] 

^ 4^2 ’ 

pL 

M3 = / exp(-^)exp[2ifc'z2]dz2 
Jo " 

1 — exp[2ifc'L] exp[—L/P\ 

1 - 2 ik'J ■ 


(A 8 ) 


(A9) 


(AlO) 
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